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Abstract: The stability problem for singular neutral functional differential systems is discussed in this 
paper. By applying a new V-functional method and the stability of the difference operator, 
some asymptotic stability criteria are derived for neutral singular differential systems with 
multiple time-varying delays. The criteria are described as matrix equations or matrix 
inequalities, which are flexible and efficient in numerical experiments. Some examples are 
given to illustrate the results. 
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1 Introduction 


In this paper, we consider neutral singular differential systems with multiple time-varying 
delays as follows 


m m 
Ei(t) ~ Ca(t — 7) = Aa(t) + > Bix (t — ri(t)) + XO Fi(2(t — n(t))) + molet), 0) 
i=1 i=l 
where t > 0, E, C, A, B; € R"*", rank(E) =r <n. The nonlinear part of system Fy € R” 
is differential, and satisfies F;(0) = 0, j = 0,1,---,m. 7;(t) is differential and 0 < 7,(t) < 
T, T(t) < d; < 1 and to = Si inf (t — 7i(t)). Suppose F;(-) (j = 0,1,--- ,m) is the higher 
ZYU,IH1,°°° m 
order term in (-), that is 


m IEE) 


=0, j=0,1,---,m. 2 
IIo Jz] f (2) 





Recently, the stability analysis of time-delay systems receives much attention of researchers. 
Many excellent results for stability of the systems with delay have been obtained!!"1) 

Singular differential equations with delay play important roles in mathematical model- 
ing of real-life problems arising in a wide range of applications, for example, multibody me- 
chanics, prescribed path control, electrical design, chemically reacting systems, biology and 
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biomedicine!!!12], But few studies on stability of singular differential equations with delay 
have been conducted so far. 

The purpose of this paper is to establish the asymptotic stability criteria for system (1). 
The complex nature of singular neutral functional differential systems causes many difficulties 
in the analytical treatment of such systems. We introduce the stability of the operator D(a.) = 
Ex(t)~Ca(t—7) and investigate V-functional method of the general singular differential delay 
system to derive the asymptotic stability criteria. The stability results have been described as 
matrix equations or matrix inequalities which are computationally flexible and efficient. 


2 Stability of the difference operator 


In this section, we introduce the stability of the difference operator. 
To simplify the discussion, the operator D : C([—7, 0], R”) — R” is defined by 


D(a) = Ex(t) — Ca(t — 7). (3) 


Thus system (2) can be written as follows 
d 
qP) = Arlt) )+ Do Bal- (0) +D Rll z(t — ri(t))) + Fo(x(t)). (4) 


Definition 2.10 The operator D is said to be stable if the zero solution of the homoge- 
neous difference equation 


D(a:)=0, t>0, to=ye {be C([-7,0],R"): Dy =0} (5) 


is uniformly asymptotically stable. 
Lemma 2.1!" The operator D : C([~r,0],R”) — R” is defined to be 


D(a) = x(t) — Cx(t — 7), 


if ||C|| < 1, then the operator D is stable, where || - || is any matrix norm. 
Since rank(E) = r < n, then it is easy to find there exist nonsingular constant matrices Pı 
and Qı, such that 
1 0 Ci Dı 
P EQ = » PCQ=|]_ _ 
0 0 D2 C2 
When detC, # 0, there exist nonsingular matrices P and Q such that 


E, 0 C 0 
PEQ = » PCQ= (6) 
0 0 0 Cy 


where C2 = Ĝi — D,C, "Da, FE, is an r x r nonsingular matrix, Cı and C are r x r and 
(n —r) x (n—r) constant matrices. 
Lemma 2.25 The operator D is stable if E7'Cy|| < 1 and detC, Æ 0. 
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3 V-functional method for singular FDEs 


In this section, we will introduce the V-functional method for general singular functional 
differential equations. 


Consider the general singular functional differential equation (FDE) as follows 


where f(t,y) : Rx C = R”, Č = C([-7,0],R”), and f(0,0) = 0, z:(80) = z(t +8), 0 € 
[-7,0], E € R”*” and rank(E) < n. For ọ € C, we define 


lel = sup lp). (8) 
0€[—7,0] 

In this paper, we always suppose that the conditions to the existence of the solutions for 
system (7) are satisfied. 

Lemma 3.1 For system (7), if there exists a definite V-functional V(y), y € C, such that 
V(y) is a semi-definitive functional with the symbol that opposite to the symbol of V(y), or 
V(y) identically equals to zero, then the zero solution of system (7) is stable. 

Lemma 3.2 For systems (7), if there exists a definite V-functional V(y), p € C, such 
that V(v) is a definite functional with the symbol opposite to the symbol of V(y), then the 
zero solution of system (7) is asymptotically stable. 

Lemma 3.3 For any z, y € R”, y > 0, inequality 2x7y < yal x + iyTy holds. 


4 Asymptotic stability criteria for singular NFDE 


In this section, some asymptotic stability criteria for system (1) are derived by applying the 
V-functional method and stability of difference operator. 

The necessary condition for the stability of system (1) is that the operator D be stable!!9, 
The following theorem governs the system (1). 

Theorem 4.1 Suppose the operator D is stable and B; (i = 1,2,--- ,m) are nonsingular. 
The zero solution of system (1) is asymptotically stable, if there exists a symmetric and positive- 
definite matrix P such that the following matrix equation 

= 1 
ATPE + ETPA (cBT PB: -pn E" PE) + ETQE =0 9 
+ +) (eB; + a) + ETQ (9) 


i=1 


holds, where € is a given positive number and Q is a given positive-definite matrix. 


Proof Construct V-functional as follows 
Vit, Zz) = Vi (t, x) + V(t, Xt); 
Vi (t, x+) = (Ex(t))? PEz(t), (10) 
Va(t, £) = € 2 ae xT (s)B? PB;x(s)ds. 
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Since B; are nonsingular and P > 0, we know BT PB; > 0. Thus V(t, x+) is positive-definite. 
A straight-forward computation gives the time derivative of the V-functional as 


Vi (t, 4) = 27 (t)(AT PE + ET PA)a(t) + 2x27 (t)ET PFo(t, 2(t)) 


Haar (t)E? PB,x(t — 7:(t)) ) +25 a" )ET PF,(t, x(t — 7;(t))), 


i=l 


no(t, 24) = =e drat (t)B? PB,x(t) epee — 7;(t)) B7 PBa(t — rit) (1 — F(t). 


t=1 
Since P is a positive-definite symmetric matrix, it has a decomposition P = UTU. From Lemma 
3.3, we get 
2a" (t)E” PByx(t — ri(t)) = 2x7 (t)ETUTU Biz (t — 7i(t)) 


1 


1 


F yet oa) a” (t)E* PEx(t) + ae(1 — di)x” (t — ri(t)) By PBix(t — 7:(t)), 


27 (t)ET PFo(t, x(t)) < 3" QE" P*Ex(t) + BFE (t,2(t)) Fo(t, 2(t)), 


2a" (t)E” PF,(t, c(t — 7i(t))) < WEY (t, a(t — ri(t))) Fi (t, e(t — ri(t))) + <2" (NET P? Ealt), 


i 


where a, 8, yi (i = 1,2,--+ ,m) are positive numbers. From above, it can be obtained 
Vem) < 27(t)[- EDE ay (G- ETP E+ G Ee L EP Eļ|z (t) 
+ Da = 1)e(1 = dija? (t = Ti(t)) B7 PByx(t = ri(t)) 
i=1 


+30 nFF (t, 2(t — 7i(t))) F(t, e(t — T:(t))) + BFS (t, x(t) Fo(t, x(t). 


i=1 
We can choose candidate positive number a(< 1), 8, yi (i = 1,2,- -- ,m) with 1—a sufficient 
small and 8, y; sufficiently - such that 


=-0+ > Gray a 1)ETPE + G +o a (11) 
Since 0 < a@ < 1 and B; (i =1,2,--- ,m) are nonsingular, we obtain 
N; = (a — 1)e(1 —d;)B] PB; <0, i=1,2,---,m. (12) 
Then we have 
V(t, z6) < 27 (t)Noz(t) Soi T(t — nE)Niz(t— rilt)) 
i=1 


HYO uF? (z(t — ri(t))) F: (z(t — r: (t))) + BEY (2(t)) Fo (z(t). 
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From (11) and (12), we can choose o > 0, such that 
N;+oI <0, j=0,1,2,---,m. (13) 


From (2), there exists 6 > 0, when |jæ(t)|| < 6, t > to, the following inequalities hold simulta- 


neously 


Poe sË < Blo, 


|R rt- OI? < Zle ETONE (14) 
Then we obtain 


m 
V(t, 4) < 27(t)(No t+ ol)a(t) + X aT (t - r(t) (Ni + ol)a(t — ri(t)). 

i=1 
From (13) we know V(t, x+) is negative definite. Since the operator D is stable, according to 
Lemma 3.2, we know that the zero solution of system (1) is asymptotically stable and this 
completes the proof. 

Theorem 4.2 Suppose the operator D is stable, the zero solution of system (1) is asymp- 

totically stable, if there exists symmetric and positive-definite matrices P and Q such that the 
following matrix equation and inequality hold simultaneously 





ATPE + ETPA = -(m + 1)ETQE, as 
a let PB PE no 
i=1 i 


Proof Similarly to the proof of Theorem 4.1, the result follows from the V-functional 
method and stability of difference operator. 
If m = 1 and F; = 0 (j = 0, 1) hold, system (1) reduces to 


E(t) — Ci(t — T) = Ax(t) + Bz (t — r(t)). (17) 
Furthermore, if t — r(t) = gt (0 < q < 1 is constant) holds, system (17) is written as 
Ea(t)~ Ca(t — T) = Ax(t) + Ba(qt). (18) 


From Theorem 4.1 and Theorem 4.2, we have: 
Corollary 4.1 Suppose the operator D is stable, the zero solution of system (17) is asymp- 
totically stable, if there exists matrices P > 0, Q > 0, and a positive parameter € such that 


one of the following two conditions is satisfied 


1 
e(1 — d) 





(Hı): A?TPE+E™PA+eBTPB+ ETPE + ETQE =0, 


ATPE + ETPA = —2ET QE, 


(H2) 
||ETPB|| <v 1- d Xmin(Q). 
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Corollary 4.2 Suppose the operator D is stable, the zero solution of system (17) is asymp- 
totically stable, if there exists matrices P > 0, Q > 0, and a positive parameter € such that 


one of the following two conditions is satisfied: 


1 
(H3): ATPE+ETPA+eBTPB + gE PE + ETQE =0, 


ATPE + ETPA = —2ET QE, 


(H4) 
||ETPB|| < Va Amin(Q). 


5 Illustrative examples 


In this section, we give some examples to illustrate the derived stability results. 
Example 5.1 Consider the following second-order singular differential system with time- 

varying delay 
i(t) — $4 (t — r(t)) = -r1 (t) + 4x1 (t — r(t)) + x7 (t)ao(t) + xf (t), 


—ao(t — r(t)) = z2(t) + af (t). 


(19) 


First, since Ey = 1, C1 = $, Cz = 1, it is easy to know that ||E7*C1|| < 1 and detC2 # 0. 
From Lemma 2.2, we know the operator D is stable. Secondly, since T(t) = 0.4t, +(t) = 0.4, the 
matrices P and Q and the positive number d in the condition (H2) of Corollary 1 are chosen 


as follows 


Then we obtain 
ATPE + ETPA = —2ETQE, ||E™PBI| < V1 —dXmin(Q). 


Thus, the condition (H2) of Corollary 1 is satisfied. The zero solution of system (19) is asymp- 


totically stable. 
Example 5.2 Consider the following third-order singular differential system with time- 
varying delay 





ialt) — $a, (t — r(t)) = -x1 (t) + o(t) + xı (t — r(t)) + x2 (t — r(t)) + sin (z1 (t — 7(t))), 
ialt) — a(t — T(t)) = —221(t) — 4x2 (t) + x2 (t — r(t)) + sin [ai (t — 7(t)) + z2(t — 7(t))], (20) 


z3(t) — £3 (t — T(t)) = z1 (t) + z2(t) — z2 (t — T(t)) + sin (z1 (t)) + sin (x3(t — 7(t))). 


First, since 
1 0 0 

E = $ Ci = , C2 = 1, 

0 1 


mE 


1 
2 


it is easy to know that |E] Ci|| < 1 and detCz 4 0. From Lemma 2.2, we know that the 


operator D is stable. Secondly, since 
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the matrices P, Q and the positive number £, d in the condition (H1) of Corollary 1 are chosen 


as follows 
1 0 0 1 0 0 
1 
d=>5, e€=1, P=/]0 10], Q=]0 3 0] >0. 
0 0 1 0 0 2 


Then we obtain 


ATPE +ETPA+ eB" PB + 27 B"PE+ E"QE=0. 


Thus, the condition (H1) of Corollary 1 is satisfied. The zero solution of system (20) is asymp- 
totically stable. 
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